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SUMMARY 


The aerodynamic damping of a flexibly mounted aeroelastic model 
with a blunted conical nose and a cylindrical afterbody was measured at 
Mach numbers from O. 7 O to 1.20 at several levels of dynamic pressure 
and two weight conditions and at Mach numbers from I .76 to 2.87 at one 
weight condition. The first two free-free flexible modes of vibration 
were investigated. Also investigated at Mach numbers from 0.9 to 1.2 
was the aerodynamic damping in the first free-free modes of a model 
which had a "hammerhead” nose (the base diameter of the blunted cone 
was greater than the diameter of the afterbody which necessitated a 
reflex angle downstream from the cone base). 

Two basically different methods, the "electrical power-input" and 
the "decaying oscillations" methods were used to determine the damping 
and frequencies. The experimentally determined values are compared with 
some applicable theories. The results of the investigation indicate 
that the aerodynamic damping in the elastic modes of vibration was small 
for all configurations tested. The maximum aerodynamic damping measured 
in the first mode was on the order of 60 percent of the structural 
damping. The aerodynamic damping was found to be even less for vibration 
modes higher than the first. Reduced-frequency effects were found to be 
negligible for the range investigated. Agreement of calculated aerody- 
namic damping derivatives with the experimental results was not good. 
Generally, the experimentally determined derivatives were larger than 
those predicted by the various theories used. The Bond-Packard theory 
appeared to give the best agreement for the first free-free vibration 
mode but gave the worst agreement for the second mode. Measurements 
made on the configuration that had a hammerhead nose indicated small 
negative aerodynamic damping in the Mach number range from 0.95 to 1.00. 
Aerodynamic stiffness effects were found to be small and within the 
experimental scatter. (Wind-on frequency determination was accurate 
only to approximately 1 percent.) 
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INTRODUCTION 


The problem of the determination of oscillatory aerodynamic deriva- 
tives of aircraft and/or aircraft components has been met to some extent 
by the use of theoretical estimates and a number of derivative measure- 
ments. However, comparatively little work has been done on this problem 
as related to slender bodies of revolution representative of guided 
missiles and space launch vehicles, particularly in the experimental 
field. Launch vehicles, especially, are structures of maximum design 
efficiency with a minimum of weight available for supporting members. 

The resulting structures are relatively flexible and could be subject to 
aeroelastic instabilities particularly in the cases where lifting pay- 
loads are involved. The need for suitable theoretical approaches and 
experimental data on the oscillatory aerodynamic derivatives of such 
slender bodies of revolution is thus becoming more urgent. The lack of 
experimental data, both as an indication of the magnitudes of the deriva- 
tives involved for flexible slender bodies of revolution and as a basis 
for evaluating the accuracy of theoretical predictions, is therefore of 
increasing concern. 

It is for this reason that the present investigation was undertaken 
to attempt to measure experimentally the aerodynamic damping and stiff- 
ness of an aeroelastic model having scaled physical properties approxi- 
mating those of a launch vehicle. The purpose of the investigation was 
twofold: first, to attempt to apply some known general techniques of 

measurement of oscillatory aerodynamic derivatives to the rather special- 
ized problem at hand; and second, to present the results of the inves- 
tigation as an indication of the relative accuracy of several theoretical 
approaches that might be used. 

There are, of course, many methods that may be employed for the 
measurement of aerodynamic derivatives in wind tunnels. Reference 1 
presents a rather complete discussion of various methods that have been 
used in the past and contains an extensive list of references on the 
subject. For the present investigation, it was decided that, because 
of the expected small values of the derivatives with the attendant dif- 
ficulty of accurate measurement, two basic techniques would be used 
concurrently. One technique chosen was the "power-input” method where 
the model is excited at resonance with an electromagnetic shaker. The 
other technique chosen was the well-known "decaying oscillation" method. 
In both techniques the aerodynamic stiffness is indicated by the change 
in frequency between the wind-on and the wind-off conditions. 

Tests were conducted in the Langley 8 -foot transonic pressure tun- 
nel over the Mach number range from 0.70 to 1.20 and in the Langley 
Unitary Plan wind tunnel over the Mach number range from I .76 to 2 . 87 . 

The basic model was a slender cone cylinder designed to approximate the 
geometric and dynamic characteristics of a representative launch vehicle. 
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Tests were also conducted on a model having natural frequencies approxi- 
mately twice those of the basic scaled model and on a model where the 
conical nose was replaced with a "hammerhead" type nose. 

Several aerodynamic theories (refs. 2 to 7 ) were used to calculate 
predicted aerodynamic damping coefficients. The results are compared 
with the experimental results. 


SYMBOLS 


Ac 

distribution function of damping portion 

of aerodynamic . 

Al 

distribution function of inertia portion 

of aerodynamic 

a k 

distribution function of stiffness portion of aerodynami 
loading 

a(t),b(t) 

tip amplitudes 


B,C 

generalized damping 


C 

damping per unit length 


^cr 

critical value of damping 


Ch 

damping derivative (c^ = 2pk-^| 

\ ^cr/ 


O 

9 

section steady-state normal-force slope 


EI 

flexural stiffness 


F(t) 

sinusoidal driving force 


F 

generalized aerodynamic force 


Gi,G2 

functions of and 0 


h 0 

amplitude of decaying vibration 


h P 

amplitude of decaying vibration after p 

cycles 

h(x),H(x) 

mode shapes based on unit tip deflection 

I 

current 
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K 

k l' k 2 

L 

M 

Ko 

m(x) 

m,n 

P 

P 

q. 


<ii><i2 

Qh 

R 

R(x) 

r 

t 


T 

V 

W(x,t) 


x,y 

y(x,t) 

Xi,X 2 

x 5 


generalized stiffness 
spring constants 

length of model 
generalized mass 
free- stream Mach number 
mass per unit length 
integers 

electrical power 

number of cycles between amplitudes h 0 and h p 
dynamic pres sure ; 

proportionality constants 

generalized force associated with aerodynamic damping 

electrical resistance 

function defining model geometry 

radius of cylindrical portion of model 

time 

arbitrary period of time 
free-stream velocity 
aerodynamic running shear load 
Cartesian coordinates 
deflection shapes 

longitudinal locations of springs k x and k 2 , respectively 

longitudinal location of applied sinusoidal force, F(t) 
angle of attack 


a 
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5b virtual displacement 

5W ^ virtual work done by aerodynamic forces 


5 ( x— x n ) 

7 

0 

P 

P 

A ,n,<u 


Dirac function 
ratio of specific heats 
cone semivertex angle 
test medium density 

/ M. \ 

mass ratio p = — - — 

\ pjcr^L/ 
circular frequencies 


Subscripts: 


m 

mth 

natural 

vibration 

mode, 

m = 1, 

2, 

3, • • 

n 

nth 

natural 

vibration 

mode, 

n = 1; 

2, 

3, • . 


A aerodynamic 

r random quantities 

s structural 

Dots over symbols indicate derivatives with respect to time. 


APPARATUS AND TEST PROCEDURE 


Model and Support System 

Construction . - The construction of the model and sting used in the 
tests is illustrated in figures 1, 2, and 3« The model was built in 
three sections to make insertion of instrumentation and ballast weights 
easier and to permit easier assembly on the sting. Basically, each 
section consisted of a fiber-glass shell 0.03 inch thick bonded to 
aluminum rings which served to stiffen the model in the hoop plane . The 
rings also were used as mounts for the accelerometers, shaker moving 
coils, the flex springs which supported the model on the sting, and the 
lead weights which were used to lower the natural vibration frequencies 
to scaled values approximating those that might be considered typical 
of launch vehicles. Model dimensions are presented in figure 2. 
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The model was supported on the sting by the flex springs shown in 
figure 5 . These springs were located approximately at the first free- 
free node points on the model in order to minimize their influence on 
the free-free modes. 

The model was excited in its natural free-free modes of vibration 
by the electromagnetic shaker built into the sting. Also incorporated 
in the sting was a provision for cooling the shaker by means of circu- 
lating cooling air through the sting and blowing it back on the field 
coils. The sting also carried a set of pneumatic brakes which were 
used to restrain the model during starting and stopping transient flow. 
These features are shown in figure l(b) . The basic model mounted in 
the Langley Unitary Plan wind tunnel is shown in figure 4. 

Model physical properties .- Three model configurations were tested. 
Configurations 1 and 2 differed only in the mass, the model geometry 
being the same. Configuration however, consisted of the basic model 
with a large hammerhead nose which was fitted over the basic shell in 
such a manner that it touched the basic shell only at the nose and 
shoulder. (See fig. 2.) The experimentally determined frequencies, 
mode shapes, and damping ratios of the three configurations are pre- 
sented in table I and figure 5* Also presented in figure 5 are the 
calculated mode shapes for a uniform free-free beam. As is shown by the 
figure the mode shapes for the models tested in this investigation are 
very similar to those of a uniform free-free beam. These mode shapes 
were determined through the use of a small velocity pick-up. The mass 
distribution, determined by weighing the individual components before 
assembly, is shown in table II. Also shown are the total weights of the 
conf igurations, the calculated generalized masses, the natural frequen- 
cies of vibration, and the values of critical damping. 


INSTRUMENTATION 


The model instrumentation consisted of an accelerometer located in 
the nose which was used to indicate the shaking amplitude of the model; 
two skin thermocouples located approximately at one-quarter and three- 
quarters of the model length; and electrical wire resistance type strain 
gages located on the flex springs. Thermocouples also were imbedded in 
the moving coils and field coils of the shaker in order to provide a 
means for monitoring the coil temperatures. An accelerometer on the end 
of the sting provided a means of determining vertical sting motion. 

A variable-frequency oscillator and power amplifier were used to 
provide alternating current to the moving coils and a direct-current 
variable voltage transformer supplied the direct current to the shaker 
field coils. A fouling light provided a means of determining that the 
moving coils were riding free of the stationary parts of the shaker. 
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The readout instrumentation consisted of the following: 

Vacuum tube voltmeter: indicated model vibration amplitude from 

accelerometer in mode nose 

Ammeter: indicated current applied to moving coils of shaker 

Ammeter: indicated current applied to field coils of shaker (held 

constant at 2.5 amps) 

Low-power-factor wattmeter: indicated power drawn by moving coils 

Oscilloscope: provided means for monitoring wave form of force 

output of shaker and of model response 

Electronic period counter: used to determine frequency of forced 

oscillations 

Potentiometer: indicated temperature of the model skin, the moving 

coils, and the field coils 

Electronic dampometer: measured the logarithmic decrement and the 

decaying frequency and was used to monitor the damping in order to 
determine any large changes due to structural failure. This electronic 
device is described in reference 8. 

Oscillograph: recorded the signals from the accelerometer in the 

model nose, the strain gages on the model supporting flex springs, and 
the fouling indicator, 

A block diagram of the instrumentation setup is shown in figure 6. 

A variable -frequency voltage was applied to the shaker moving coils by 
an oscillator and a power amplifier. A constant current of 2.5 amperes 
was supplied to the shaker field coils by a d-c power supply. The cur- 
rent through the moving coils was measured by use of a vacuum-tube 
voltmeter shunted with a 1 -ohm precision resistor. The power consumed 
by the shaker was measured by use of a low -power -fact or wattmeter. A 
supplementary circuit was employed to determine accurately the resistance 
of the moving coils which varied due to temperature changes in the coils. 
Immediately after each test point, a battery- supplied voltage was applied 
to the coils and the voltage and corresponding current noted from the 
voltmeter and ammeter, respectively. 


Wind Tunnels 

The investigation was conducted in two wind tunnels. The Langley 
8-foot transonic pressure tunnel was used for the Mach number range from 
0.8 to 1.2 and the Langley Unitary Plan wind tunnel was used for the 
Mach number range from 1 . 7 6 to 2.87. Both tunnels are of the continuous- 
flow type with variable Mach number and density control. Stagnation tem- 
perature in both tunnels was held constant at approximately 125° F . 
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Test Procedure 

The test procedure used in both tunnels was essentially the same. 
After mounting the model in the tunnel it was shaken at the resonant fre- 
quency of the desired free-free mode with no air flowing at several test- 
section pressures from atmospheric down to moderately low vacuums with 
cooling air on and off. No change in frequency or damping was observed. 
While the model was vibrating at a given amplitude, data were recorded 
from all the meters necessary for application of the "power-input" method 
of determining damping. The power to the shaker was then cut (and the 
shaker moving coils removed from the electrical circuit) and the decaying 
oscillations of the model were recorded. While the model motion was 
decaying, the logarithmic decrement and decaying frequency were measured 
by an electronic dampometer. This device provided a quick means of 
detecting any extreme changes in damping due to structural failure, loose 
screws, and so forth. After obtaining the no-wind data, the model "brake" 
was actuated and the tunnel started. When the tunnel reached the desired 
operating conditions, the model brake was released and the above procedure 
repeated. The process was continued for all the modes investigated at the 
various Mach numbers and at several levels of dynamic pressure at each 
Mach number. The Reynolds number range for the tests was from 0. 7 X 10^ 
to 2 X lO^ based on model afterbody diameter. 


Reduction of Data 

The method of calculating damping from freely decaying motions of a 
structure is well known. The ratio of total damping to critical damping 
for a given mode of vibration of a linear system is 

, _k_ log. 5a 

C cr 2rtp hp 


where 


*0 


an amplitude of vibration 

amplitude of vibration after p cycles 

number of cycles between amplitude h.Q and hp 


The aerodynamic damping ratio 


C A 


is determined by measuring the struc- 


■'cr 


tural damping ratio — 

^cr 


in still air before each test and subtracting 
it from the total damping measured during the test. 
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Relationships for determining the damping ratio from the elctrical 
measurements are developed in appendix A. The resulting expression is 

C A + C s = ^-1 2 h ( x 5) I 2 
C cr 2q2^s^ P - I^R 


(see appendix A) where 


H 

I 

M s 

P 

*1 




R 

A 


mode-shape deflection at shaker location 
current through shaker moving coils 
generalized mass of model 

electric power consumed by shaker moving coil circuit 

proportionality constant relating shaker force output to the 
current applied to moving coils 

a proportionality constant relating model tip velocity to 
voltage applied to moving coils 

electrical resistance of shaker moving coil circuitry 

natural vibration frequency at test condition 


Since the ratio of structural damping to critical damping of the model 
was found to be slightly dependent upon amplitude and temperature, it 
was necessary to construct calibration curves of structural damping 
ratios against model- vibration amplitude and skin temperature so that 
structural damping measured before the tests could be corrected to 
correspond to that associated with the vibration amplitude and skin 
temperature of the model at the test point, A representative set of 
these calibration curves is shown in figure 7- After adjusting the 
no-wind structural damping to the value associated with the wind-on test 
point, the adjusted structural damping ratio was subtracted from the 
wind-on total damping measured to get the aerodynamic damping ratio. 

For convenience in presenting the data, an averaged value of the struc- 
tural damping ratio was then added back to the aerodynamic damping ratio 
for each test condition to give the total damping ratio that is presented 
in the results. 


The ratios of aerodynamic damping to critical damping obtained in 
this manner apply only to the particular mechanical systems used in these 
tests. In order to make the damping measurements of more general value, 
a nondimensional damping force coefficient C£ can be defined in terms 
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of the usual flutter parameters, reduced frequency and mass ratio, in 
such a manner as to permit estimation of aerodynamic damping of similar 
configurations having somewhat different structural characteristics. 

As indicated in appendix B (eq. B5) 



p 

irtpVr 

2M S A 



and is defined by 

C A _ JtpVr 2 c . 
C cr 2M s A h 


or, if this expression is rewritten in terms of conventional flutter 
parameters. 



where 

k 

A 

V 

r 

P 

L 

M 

1A s 

A 


reduced frequency, 

Mo 

an effective mass ratio, — - — 

p 

nr pL 

flow velocity 

radius of cylindrical portion of model 
mass density of test medium 
overall length of configuration 

generalized mass of configuration in vibration mode of 
interest 

frequency of vibration in mode of interest 


This damping-force coefficient C£ is useful in scaling model results 

to equivalent full-scale conditions and is used in the presentation of the 

C A 

results of the investigation along with the damping ratio — — . 

C cr 
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RESULTS AND DISCUSSION 


The values of damping observed in this investigation for the various 
test conditions are shown in figures 8 to 11. Shown in figure 8(a) is 
the variation of the ratio of total damping to critical damping with the 
flow parameter pV (test- section air density multiplied by the free- 
stream velocity) for model configuration 1 oscillating in the first free- 
free mode at several Mach numbers. The parameter pV is used because the 
theories applied have all predicted a linear variation of aerodynamic 
damping with pV for a given Mach number. (See appendix B.) The dashed 
line is the averaged structural damping ratio. Shown plotted are data 
obtained from the power input method (open symbols) and from freely 
decaying oscillations (solid symbols). The freely decaying oscillations 
data is an average of that obtained from the electronic dampometer and 
from the oscillograph records. It should be mentioned that the free-decay 
method could not be used transonically because of the very uneven decays 
caused by the model response to the transonic flow conditions. Although 
there is scatter in the data, some observations can be made. First, note 
that the aerodynamic damping ratio is small, being on the order of 60 per- 
cent of the structural damping ratio at the maximum pV tested. 

Also, note that for the supersonic Mach numbers investigated in the 
Langley Unitary Plan wind tunnel, the ratio of aerodynamic damping to 
critical damping appears to vary linearly with pV whereas in general 
this is not true for the transonic Mach numbers. 

In figure 8(b) the aerodynamic derivative for the first free- 

free vibration mode of model configuration 1 is plotted against Mach num- 
ber for several of the higher values of pV. Data from both the power- 
input and free-decay methods of measurement are included for the Mach 
number range from I .76 to 2 . 87 . Within the scatter of the data (which is 
greatly magnified on the basis of the parameter C£), there is apparently 

little Mach number effect. Also shown in figure 8(b) are some aero- 
dynamic damping derivatives calculated (using experimentally determined 
mode shapes) from piston theory, momentum theory, second-order shock- 
expansion theory, second-order Van Dyke theory, and a theory by Bond and 
Packard. (See appendix B for an explanation of how these various theories 
were applied.) The second-order shock-expansion and second-order Van Dyke 
theories show the worst agreement with experiment. The trend with Mach 
number indicated by piston theory is not evident in the experimental 
results. The Bond-Packard theory agrees very well with the experimental 
results . 

The results of damping measurements in the second free-free vibra- 
tion mode of model configuration 1 are shown in figure 9* The variation 
of total damping with the flow parameter pV at several Mach numbers is 
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shown in figure 9 ( a ) • Note that, unlike the first-mode data, the struc- 
tural damping as determined from the power input and free-decay methods 
is not the same. The reason for this wide disparity has not been 
resolved and warrants further investigation. It is felt, however, that 
under forced vibration conditions the model may be responding in a mode 
that degenerates into a slightly different shape when the model is freely 
decaying. If this is the case, the change is not significant enough to 
alter the trend reflected in the aerodynamic damping shown for the two 
methods in figure 9 ( a )* Although the structural damping differs almost 
by a factor of two, the aerodynamic damping measured by the two methods 
is in fairly good agreement. In both cases, the maximum aerodynamic 
damping is less than k-0 percent of the corresponding structural damping. 

In figure 9(^0 j the damping derivative C£ for the second free-free 
mode is plotted against Mach number for several of the higher values of 
pV. Experimental values determined from the two methods are shown with 
some calculated derivatives. No Mach number effect is apparent from the 
power- input data although the free -decay-method data show an apparent 
decrease in the damping derivative with Mach number over the range that 
the method was usable. The same trend is seen to be predicted by the data 
calculated by using the piston theory. The other theories are not in good 
agreement with the experimental values . 

The results of damping measurements on the first free -free mode of 
model conf iguration 2 (which had a frequency about twice that of con- 
figuration l) is shown in figure 10. The variation of total damping 
ratio with the flow parameter pV is shown in figure 10(a) for several 
Mach numbers. The primary observation to be made here is that the maxi- 
mum aerodynamic damping is about 60 percent of the structural damping, 
which is about the same as it was for configuration 1. Also, the values 
of C-^ for the two configurations are roughly the same at comparable 

Mach numbers. The results of tests on the two configurations therefore 
indicate little, if any, reduced-frequency effect for the range 
investigated. 

In figure 10(b) the variation of the aerodynamic damping deriva- 
tive C£ with Mach number is shown for several values of pV. No Mach 

number effect is apparent in the experimental data. Also shown are 
values of obtained from quasi-steady calculations based on static 

pres sure -distribution measurements made on a rigid model of like geometry 
given in reference 9 - These calculated values are seen to be low with 
respect to the experimental damping derivatives as was the case for the 
purely theoretical approaches used on configuration 1. 

Figure 11 presents the variation of total damping with Mach number 
for model configuration 5 vibrating in the first free-free bending mode. 
The interesting point to be observed here is that the addition of the 
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11 hammerhead" nose to configuration 1 resulted in small negative aero- 
dynamic damping over a portion of the Mach number range from about O .95 
to 1.00. 

It might be mentioned that the scarcity of data points for this 
configuration relative to the other configurations tested is due to the 
fact that the model response to buffeting loads was so severe that the 
springs that supported the model on the sting kept breaking before the 
data could be obtained. 

A few remarks on the results of the test program in general are in 
order. Although an attempt was made to measure aerodynamic stiffness in 
addition to aerodynamic damping in the elastic modes, no discernible 
change in model frequency was apparent outside the normal scatter in the 
data between the wind-off and the various tunnel wind-on conditions . 
(Wind-on frequency determination was accurate only to approximately 
1 percent.) Likewise, attempts to measure aerodynamic damping in the 
third free -free mode proved to be futile because of the apparently 
extremely low values even under maximum pV conditions. Also, "hoop" 
modes appearing in the third free -free bending mode made accurate 
determination of the mode shape extremely difficult. 


CONCLUDING REMARKS 


Aerodynamic damping measurements made on a sting -mounted aero- 
elastic model with a blunted conical nose and a cylindrical afterbody 
over a Mach number range from 0.8 to 2 . 87 , and on a model with a 
"hammerhead” nose at Mach numbers from 0.9 bo 1.2 indicate that the 
aerodynamic damping in the elastic modes of vibration was small for all 
configurations tested. The maximum aerodynamic damping measured was 
approximately 60 percent of the structural damping. The aerodynamic 
damping was found to be even less for vibration modes higher than the 
first . 

Reduced-frequency effects were found to be negligible for the range 
investigated . 

Of the various theoretical methods considered, only the Bond- 
Packard theory for damping in the first elastic mode showed good agree- 
ment with the experimental results. Generally, the experimentally 
determined derivatives were larger than those predicted by the various 
theories used. 

Measurements made on the conf iguration that had a hammerhead nose 
indicated small negative aerodynamic damping in the Mach number range 
from 0.95 to 1 . 0 . 
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Aerodynamic stiffness effects were found to be small and fell within 
the experimental scatter since wind-on frequency determination was accu- 
rate only to approximately 1 percent. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 7> 1962. 
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APPENDIX A 

DEVELOPMENT OF DATA -REDUCTION RELATIONS 
FOR FORCED OSCILLATION DAMPING MEASUREMENT TECHNIQUE 


The purpose of this appendix is to present in some detail the rela- 
tionships necessary to the data reduction for the forced-oscillation 
damping-measurement technique . It should he pointed out that a somewhat 
similar presentation is given in reference 10 for a cantilever wing 
oscillating in the first bending mode. 

The system considered here is a flexible beam having running 
mass m(x) and stiffness El(x) which is attached to two soft linear 
springs kp and k 2 at longitudinal coordinates Xp and x 2 - The 

system is excited by a sinusoidal force F(t) applied perpendicular to 
the beam center line at longitudinal station x^. In addition to the 

force F(t) and the inertia and spring force loading due to the motion 
of the beam, the beam is subjected to a running aerodynamic shear 
loading W(x,t) . The system is presented schematically in sketch 1. 

The x,y coordinate system is fixed relative to the undisturbed beam. 



In the absence of the driving force F(t) and the aerodynamic 
loading W(x,t) the differential equation which defines the transverse 
bending free -vibration characteristics for the undamped system is 


S^i( x )Sfy(x i t)_ + k i y ( x , t ) &(x— xp) + k 2 y(x,t) &(x-x 2 ) = -m(x) d 2 y( x ' t} 
Gx 2 Sx 2 °t 

(Al) 



i6 


where the function 5(x— x n ) 


is a Dirac function with the property that 



S(x— x n ) dx 


y( x n^) 


(A2) 


The displacement function y(x,t) may be expanded into a series of prod- 
ucts of time and space functions . 


00 

y(x,t) = a^t) h m (x) 


m=l 


(A3) 


By choosing the normal modal functions of the system defined by equa- 
tion (Al) for the space functions h m (x) and applying a Galerkin proce- 
dure along with the appropriate orthogonality conditions, equation (Al) 
becomes a set of equations: 


^s,nftn(^) + ^s,m a m(t) - 0 (m - 1, 2, 3* • • •) (A4) 


where the generalized mass M s ^ m is f m(x) h m ^-(x) dx and is 

J 0 

the natural frequency of vibration of the mth mode of the system defined 
by equation (Al) . A viscous-type structural damping may be included by 
adding a term B m a m (t) . 


^s,m®m(^) + ®m^m(^) + t %i'^ni a m( ^ ) - 0 ( m - 2j 3; • • •) (A5) 


Now with the sinusoidal force F(t) = F Q e imt applied to the damped 
system at the point x = x^, the set of equations of motion is 

^s,m^m(^) + ®m®m(t) + f %t^s,in a ni(^) = f F(t) h^x) 5(x— x^) dx 

J 0 

= F(t) bju( x^) (m = 1, 2, 3, . . .) (A6) 
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The deflection shape assumed by the beam after all transient motion has' 
subsided is 


y(x,t) = }_ 
m=l 


F (t) hm( x 3 ) h m( x ) 
K 2 " ^ 2 ) M s,m + 1B m 


(A7) 


If the displacement function y(x,t) were plotted against <x>, there 
would be n peaks associated with oo = . Each of these peaks would 

have a shape very similar to that found for a lightly damped mass oscil- 
lator. For reasonably veil separated frequencies 




(A8) 


Associated with each of these resonance frequencies there will he a 

particular deflection shape. The shapes are given by equation (A7) when 
no = H n . The function y n (x,t) may be separated into the product of a 
time and space function 


y n ( x ,t) = ' b n( t ) H n( x ) 


(A9) 


The function H n (x) may be thought of as a coupled-mode shape and the 
frequency Q n as the corresponding coupled-mode frequency. These modes 
are coupled as a result of the damping forces and the shaker forces. 
Because the structural damping is small, the coupled- and uncoupled- 
mode shapes and frequencies are almost identical. These shapes and 
frequencies may be calculated or measured experimentally. The latter 
method was chosen for this investigation. 

In this development the motion of the system will be examined only 
at these coupled-mode resonance peaks. At the nth resonance condition 

the potential energy of the system is ^ n b n (t)^ where K n may be 
thought of as a generalized stiffness and is equal to 

dx + kpH n (xp ) 2 + k 2 H n (x 2 ) 2 
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The kinetic energy of the system is where M s n is 

generalized mass and is equal to 



m(x) H^x) dx 


The dissipation function associated with the assumption of viscous 
structural damping is ic s ^ n h n (t)2 where C s ^ n is the generalized 

damping and is equal to 



2 

c s,n H n ( x ) dx 


The virtual work associated with the external loads F(t) and 
W(x,t) acting on the system through a virtual displacement 5b n (t) H n (Xj) 
is 

6b n (t) F(t) H n (x 5 ) + 5 W a (A10) 


Before the virtual work done by the aerodynamic forces can be expressed 
more precisely than by 5W^, the nature of W(x,t) must be examined in 
some detail- This loading may be separated into three categories. 

The first is that due to the mean angle of attack a. Since this 
portion of the aerodynamic loading is steady state and produces only 
static deflections of the structure, it will not be treated further. 

The second is of a random type being produced by turbulence, sepa- 
rated flow, and so forth. These random forces will produce random 
deflections b n ^ r (t), velocities b n, r (t), and accelerations b n ^ r (t) 
of the beam. The virtual work done by the random aerodynamic loading is 



W r (x,t) H n (x) dx 


where W r (x,t) is the random aerodynamic running shear load. 
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The third type of aerodynamic loading is that due to the motion of 
the team itself. This portion of the aerodynamic loading may he sepa- 
rated into three parts: inertia, damping, and spring forces. It must 

he assumed that the aerodynamic inertia and spring forces are reasonably 
small relative to their structural counterparts, so that the mode shapes 
will not differ appreciably from those measured with zero wind. The 
virtual work done by these three aerodynamic forces is 


-5b n (t) b n (t) f Aj H n (x) dx - 5b n b n (t) 
J 0 



^(x) ax 


- 8b n (t) b n (t) 



H n( x ) 


where Aj n , A c n , and A^ n are distribution functions of the inertial, 
damping, and stiffness components of the aerodynamic loading. 

Combining these expressions the total aerodynamic virtual work is 


5W A = -^nh) 


M, 




(All) 


where 


M 


[ A,n = / A I,n H n( x ) ^ 


0 


A,n = / A c,n H n( x ) 


0 


KA,n = / A K,n H n ( x ) dx - 


0 


r,n 


0 


W r (x,t) H n ( x ) ^ 
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By applying Lagrange's dynamical equation of motion, the differential 
equation of motion of the system under the applied shaker and aerodynamic 
loading is 


( M s,n + m A, 


n 






+ ( K s,n + K A,n) + + F r,n (A12) 

v L _j 


Multiplying equation (A12) by F(t) and taking time averages yields 


(M s ,n + M^n)^) F(t) + (}A Btn + M A ^ n )b n , r (t) P(t) 


( C s,n + ^A,n)^n^^ + ( C s,n + C A,n)^n,r^ F ( t ) 

( F s,n + K^^hjt) F(t) + (k^ + K^ n )b^ r (t) P(t) 


= ^(* 3 ) F ( t ) 2 + F r,n F(t) (A13) 


where the bars indicate a time average over a period of time T. However, 
as T approaches co, the products 


M, 


;,n + ^,11)^,^^ PC*) 


C s,n + C A,n)^n,r( t ) 
K s ,n + ^njVrit) 


F r,n F (*) 


vanish since there is no correlation between the random part of the 
response and the sinusoidal forcing function and thus the time average 



of the product of these quantities is zero. Also, since the system is 
being forced at resonance, the force is in phase with b n and is 90 ° 

out of phase with b n and b n ; therefore the products 


(M s ,n +M A , n )b n (t) F(t) = ( k^ u + K A ^ n )b n (t) F(t) = 


0 (Al4) 


Consequently, 


( C s,n + * H n( x 5 ) F(t > ; 


(A15) 


or 


+ C A, n 


= ^n( x ^) F(t) ; 

b n (t) F(t) 


(A16) 


The critical damping of the system is defined as 


^cr,n 2 i(^s,n + M A,n)(^s,n + K A,n) 2 \i(^s,n + ^A,n) (^7) 


However, since M^ n is small relative to M s ^ n , the critical damping is 


^cr,n “ ^\i^s,n 


(A18) 


The frequency 7^ is the frequency of the system in the nth mode with 

the external aerodynamic loading present. The ratio of total damping to 
the critical damping is 


"A,n 


+ C 


s,n 


H 


n(x 5 ) F(t)' 


✓ \ / \ 


(A19) 
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If the sinusoidal force F(t) is applied by an electromagnetic 
shaker, it can be shown that the current passing through the shaker 
moving coil is proportional to the output force of the shaker. 


F(t) = q 1 I (A20) 

Also, the quantity b n (t) F(t) is proportional to the power consumed by 
the shaker moving coil less any power losses due to coil resistance. 

b n (t) F(t) = q 2 , n (p - I 2 R) (A21) 

Therefore, by using the relationship given by equation (A20), the expres- 
sion for the damping ratio becomes 


"A,n 


+ C 


s,n 


■'cr^n 


9]. 2 Xnfa) l2 

92,n(P - I 2 !*) K,n \) 


(A22) 


The aerodyn ami c damping ratio 


C A,n 

r 

°cr,n 


is determined by measuring the 


^ s n 

structural damping ratio — in still air before each test and sub- 

^cr, n 

tracting it as a tare from the total damping measured during the test. 
It should be pointed out that external forces on the model can produce 
coil velocities which induce correlated voltages and currents in the 
electrical system. These currents are kept small, however, by driving 
the shaker coil with a voltage generator that has an impedance which is 
high with respect to that of the moving coil. 
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APPENDIX B 

CALCULATED AERODYNAMIC DAMPING COEFFICIENTS 


The theoretical model used here is a slender cone cylinder oscil- 
lating harmonically about a mean angle of attack of 0° in its nth 
transverse bending mode at frequency 7^ and with tip amplitude b n (t) 

and normalized mode shape H^x) . The theoretical model differs from 

the experimental one in that the conical nose section is assumed to be 
sharp. Since only the aerodynamic damping is desired, only those forces 
generated by the vertical velocity b n (t) H n (x) are considered. 

The generalized force associated with the aerodynamic damping is 

Si, n -*n(t) / A H^x) dx (Bl) 

’ J 0 

where A c is the distribution of the aerodynamic damping forces per unit 

length. The subscript n will be omitted hereafter since reference is 
made only to a system oscillating in a single mode . 

Using a quasi-steady approach the distribution function A c may be 
expressed as 


Ac 


2qrtR(x) C K ^ a ^ 


(B2) 


where R(x) is the function which describes the shape of the body, ^N,a 
is the section steady-state normal-force coefficient slope, and is 
a normalized angle of attack. Equation (Bl) may now be written 




C N?a R(x) H 2 (x) dx 


b(t) 


(B3) 


where the term in the brackets is C^. 



2k 


C A = rt P v f c N, a R ( x ) h2 ( x ) (B4) 

J o > 

or , in nondimensional form, 



p 

rtpVr 

2M S A 



(B5) 


The theories of references 2 to 4 are readily adaptable to the quasi - 
steady approach. In general, the equation for the aerodynamic damping 
coefficient (b 4) will have to be integrated in two parts - one part for 
the conical nose portion and one for the cylindrical afterbody. The 
expressions for the section normal-force coefficient as predicted by the 
methods presented in references 2 to 4 are presented in the following 
table : 


Theory 

Reference 

c N,a for the - 

Conical nose 

Afterbody 

Second-order 
Van Dyke 

2 

2 tan 8 ^1 - tan^8 

to l°ge( , ^ W + l] 

\JVL 2 - 1 tan 0/ 


0 

Second -order 
piston 

5 

^ + ( 7 2 )tan 9 + j^ 7 ^ tan 2 Q 

1 

Moo 

Second-order 

shock expansion 

k 

Exact solution for cone (ref. 5 ) 

, -g 2 Z. 


The functions Gp and G 2 are functions of and 8 alone and X 

is a distance downstream from the cone-cylinder intersection. (See 
appendix C of ref, 4.) 

Another theory which lends itself readily to aerodynamic damping 
calculations is the linear small -aspect -ratio unsteady-flow theory given 
in reference 6. This theory gives for A c 


A c = jtpV 


2R(x) 2 -^H(x) + 2R(x) -^B(x) H(x) 
dx dx 


(b6) 
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The aerodynamic damping coefficient is 


C A = rtpV 



^[r 2 (x) H(x) 2 ]cbc 


(B7) 


Performing the integration in equation (B 7 ) yields a very concise 
expression for the aerodynamic damping 


C A = jrpVr 2 H(L ) 2 


(b 8 ) 


where H(L) is the value of the normalized deflection shape at the rear 
of the model, and r is the radius of the cylindrical portion of the 
model . 

In reference 7 Bond and Packard present an analysis based on a 
slender -body approximation, the solution of which is obtained in a defi- 
nite integral that has not been found to be integrable in finite form. 
Evaluation of the integral has been accomplished by numerical methods, 
and the computation forms have been programed for the IBM 70b data 
processing machine for certain cases. 

The analysis of reference 7 has been used by the Ames Research 
Center to calculate values of the damping force coefficient C£ for the 

first two flexible modes of model configuration 1 at supersonic Mach 
numbers . 

The results of those calculations are presented in figures 8 and 9 
along with the results of the other theoretical approaches used. 



26 


REFERENCES 


1. Bratt, J . B . : Wind Tunnel Techniques for the Measurement of Oscilla- 

tory Derivatives. British A.R.C. 22,146, Aug. 19 > i 960 . 

2. Van Dyke, Milton D.: A Study of Second-Order Supersonic Flow Theory. 

NACA Rep. 108l, 1952. (Supersedes NACA TN 2200.) 

5 . Miles, J. W., and Young, Dana: Generalized Missile Dynamics Analysis. 

Ill - Aerodynamics. GM-TR-OI 65 -OO 36 O, Space Tech. Labs., The Ramo- 
Wooldridge Corp., Apr. "J > 1958. 

4. Syvertson, Clarence A., and Dennis, David H.: A Second-Order Shock 

Expansion Method Applicable to Bodies of Revolution Near Zero Lift. 
NACA Rep. 1328 , 1957- (Supersedes NACA TN 5527-) 

5. Stone, A. H.: On Supersonic Flow Past a Slightly Yawing Cone. Jour. 

Math, and Phys., vol. XXVII, no. 1, Apr. 1948, pp. 67 - 81 . 

6 . Garrick, I. E.: Some Research on High Speed Flutter. Third Anglo- 

American Aero. Conf., Sept. 4-7, 1951 (Brighton, England). R.A.S., 
1952 , pp. 419-446J. 

7. Bond, Reuben, and Packard, Barbara B.; Unsteady Aerodynamic Forces 

on a Slender Body of Revolution in Supersonic Flow. NASA TN D-859 j> 
1961 . 

8 . Olsson, Carl Olof, and Orlik-Riickemann, Kazimierz: An Electronic 

Apparatus for Automatic Recording of the Logarithmic Decrement and 
Frequency of Oscillations in the Audio and Subaudio Frequency Range. 
Rep. No. 52, Aero. Res. Inst, of Sweden (Stockholm), 195^ • 

9- Coe, Charles F.: Steady and Fluctuating Pressures at Transonic Speeds 

on Two Space-Vehicle Payload Shapes. NASA TM X-503; 1961 . 

10. Rainey, A. Gerald: Measurement of Aerodynamic Forces for Various Mean 

Angles of Attack on an Airfoil Oscillating in Pitch and on Two 
Finite -Span Wings Oscillating in Bending With Emphasis on Damping 
in the Stall. NACA Rep. 1305, 1957* (Supersedes NACA TN 3643.) 


o o 


27 


TABLE I.- EXPERIMENTALLY DETERMINED NATURAL MODE SHAPES 
OF MODEL CONFIGURATIONS TESTED (NORMALIZED ON NOSE 
OF CONFIGURATIONS 1 AND 2) 


x/L 


.019 

.10 

.20 

• 25 
.50 

• 35 
.4o 
•50 
.60 
.65 
.70 

• 75 
.80 

.85 

•90 

1.00 


Normalized deflections 

Configuration 1 

Configuration 2 } 

Configuration 3* 



mode 1 

mode 1 

Mode 1 

Mode 2 







1.091 

1.000 

1.000 

1.000 

1.000 

.620 

.462 

.564 

.577 

.235 

-.095 

.125 

.155 

.050 

-.506 

-.095 

-.055 

-.125 

-.451 

-.303 

-.231 

-.273 

-.494 

-.472 

-.368 

-.385 

-.417 

-.596 

-.461 

-.478 

-.073 

-.712 

-.538 

-.4io 

• 335 

-.644 

-.447 

-.316 

.450 

-.553 

-.338 

-.197 

.476 

-.417 

-.184 

-.038 

.588 

-.239 

0 

.133 

.225 

-.018 

.195 

.520 

.002 

.209 

.398 

.517 

-.265 

.458 

.606 

.930 

-.950 

.895 

1.049 






(a) General viev of model and sting assembly. L-61-1191.1 



(t>) Detailed view of sting in vicinity of shaker. L-61-1192.1 


Figure 1.- Photograph of basic model and sting. 
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Confirmations 1 and 2 




Longitudinal section through nose cap of configuration 3 


Figure 2.- Geometry of configurations tested. All dimensions are in 

inches except as noted. 



Figure 3 .- Photograph of flex springs used to support model on sting. L-6I-II93 





l- 6 i- 96 i 

Figure 4.- Photograph of basic model mounted in Langley Unitary Plan wind tunnel. 


Measured 
Uniform beam 


(a) Configuration 1; first mode; average C s /C cr , 0.0084; 
model frequency, 85 cps; uniform beam frequency, 84 cps 


Measured 
Uniform beam 


Cb) Configuration 1; second mode; average 0 S /C cr , 0*011; 
model frequency, 223 c P s ; uniform beam frequency, 213 cps. 

Figure 5*- Comparison of measured mode shapes with those of uniform 

free-free beam. 






o 


Measured 
Uniform beam 



x/L 


(c) Configuration 2; first mode; average C s /C cr , 0.0094; 
model frequency, 159 cps; uniform team frequency, 171 cps 



x/L 


(d) Configuration 5j first mode; average C s /C cr , 0.0111; 
model frequency, 79 cps; uniform team frequency, 8l cps 


Figure 5*- Concluded. 



Thermocouples Potentiometer 
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Figure 6.- Block diagram of model, shaker, and readout instrumentation. 















Tip amplitude, in 

(a) Variation of structural damping coefficient with model 
amplitude for constant model skin temperature . 
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Model skin temper at lire, °F 


(b) Variation of structural damping coefficient with model 
skin temperature for constant model amplitude* 

Figure 7*- Representative structural damping calibration curves 
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3.0 \- 


O Power method 
• Free .decay 

Piston theory, Ref. 3 

Momentum, Ref. 6 

2d order shock expansion. Ref. 

2d order Van Dyke, Ref. 2 

— □ Bond-Packard, Ref. 7 
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(b) Variation of aerodynamic damping derivative with Mach number. 

C A c . _1_ 

C cr ' h 2pk* 


Figure 8.- Concluded. 




P V, lb-sec/ft 3 


(a) Variation of total damping ratio with flow parameter. 

Figure 9 • “ Damping measurements on model configuration 1 vibrating in 

the second free-free mode. 
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(b) Variation of aerodynamic damping derivative with Mach number. 

C A „ i 

C cr h 2pk 


Figure 9 *- Concluded. 
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(a) Variation of total damping with flow parameter. Power method. 


3.0 - 

O rower method 

• Calculated from experimentally 

determined static pressure distributions, Ref. 9 

2.0 ~ Momentum theory. Ref* 6 
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(b) Variation of aerodynamic damping derivative with Mach number. 

Figure 10.- Results of damping measurements on model configuration 2 
vibrating in the first free-free mode. 
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NASA -Langley, 1962 L— 20^^ 


Figure 11.- Variation of total damping with Mach number for model 
configuration 3 vibrating in first free-free mode. 



